The chemical potential for the inhomogeneous electron liquid in terms of its kinetic 
and potential parts with special consideration of the surface potential step and 

BCS-BEC crossover 



K. Morawetz* 1 '^. H. March 3 ' 4 * 5 and R. H. Squire 6 
Institute of Physics, Chemnitz University of Technology, 09107 Chemnitz, Germany 



Max-Planck-Institute for the Physics of Complex Systems, Nothnitzer Str. 

3 Department of Physics, University of Antwerp, 
4 Oxford University, Oxford, England 



01187 Dresden, Germany 



Abdus Salam International Centre for Theoretical Physics, Trieste, Italy and 
6 Department of Chemistry, West Virginia University, Montgomery, WV25136, USA 

(Dated: February 5, 2008) 

The chemical potential fj, of a many-body system is valuable since it carries fingerprints of phase 
changes. Here, we summarize results for /i for a three-dimensional electron liquid in terms of average 
kinetic and potential energies per particle. The difference between p and the energy per particle 
is found to be exactly the electrostatic potential step at the surface. We also present calculations 
for an integrable one-dimensional many-body system with delta function interactions, exhibiting a 
BCS-BEC crossover. It is shown that in the BCS regime the chemical potential can be expressed 
solely in terms of the ground-state energy per particle. A brief discussion is also included of the 
strong coupling BEC limit. 
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The dependence of the chemical potential ll at T = 
on a physical observable such as particle number den- 
sity is considered for two quantum liquids. As a starting 
point, the homogeneous electron liquid in three dimen- 
sions is discussed by combining the quantum-mechanical 
virial theorem with results of Hugenholtz and Van Hove 
and the Budd- Vannimenus theorem. As the electron den- 
sity is lowered, ll is related to surface properties. 

Secondly, consideration is given to a recent one- 
dimensional model possessing a crossover from a 
Bardeen-Cooper-Schrieffer (BCS) state to a Bose- 
Einstein condensate (BEC). This crossover occurs at 
fj, = 0, changing from positive values of [i in the BCS 
phase through zero into the BEC phase at negative Li. Us- 
ing the thermodynamic relations like fi = dE/dN where 
E is the ground-state energy and N is the number of par- 
ticles, E and ll are readily related in the one-dimensional 
model. 

One of us [1] in early work combined a many-body 
result of Hugenholtz and Van Hove [2| with the virial 
theorem [j| E| for a homogeneous correlated electron liq- 
uid. This led to an expression for the chemical potential 
\x as a linear combination of the average kinetic energy 
(T), and the corresponding average potential energy (V) 
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where N denotes the number of electrons. Of course 
(T) and (V) depend on the value of the uniform electron 
density n or equivalently on the mean inter-electronic 
spacing r s 
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One knows from free Fermi gas theory that as r s — ► 0, 
i.e. in the extreme high-density limit, the kinetic en- 
ergy per electron is 3/5 of the Fermi energy, which in 
this limit is the chemical potential. Thus one can ob- 
tain insight as r s — > into the exact result (|TJ) for /i in 
the correlated uniform electron liquid of density n. The 
quantum-mechanical averages of the kinetic and poten- 
tial energies in ([1]) respectively depend on n or by J5J on 
r s . As r s — > 0, the kinetic energy per particle (T)/N is 
proportional to r\T 2 and the potential energy per parti- 
cle (V) /N to rj 1 and hence the kinetic energy dominates 
the potential term, as required in the free Fermi gas limit. 
Consequently, as r s is increased, the chemical potential 
ll will pass from a large positive value at high density, 
through zero, to a negative value at sufficiently low den- 
sity. 

It is now interesting to note that the exact relation ((T|) 
can be rewritten into another form revealing the link to 
surface properties. Therefore we repeat the steps leading 
to (TT]) but with the explicit inclusion of the surface as- 
suming for simplicity that we have a flat surface and the 
electron liquid in the half space x < 0. Then the Budd- 
Vannimenus (BV) theorem [5( links the surface potential 
step of the extended electron system directly to the bulk 
energy E^/N per particle 



A$ = 0(0) 
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The electrons leak out of the surface which leads to a 
depletion of electron density inside the electron system 
such that the surface dipole develops that results in the 
electrostatic potential step ((3]). This BV theorem for 
special geometries can be found in jf| and has played an 
important role to clarify discrepancies in observations of 
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thermodynamic corrections in the Bernoulli potential at 
superconducting surfaces @, 3 ■ 

Compared to the original form Q, the virial theorem 
now contains an additional surface term Q 
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where A is the surface area and a the surface energy. 
Distinguishing the surface and bulk parts in the kinetic, 
T = T b + T a A/N, and potential energies, V = V b + 
V a A/N one obtains the virial theorems separately for 
bulk and surface properties, respectively [9( 



2(T b ) + (V b ) = -Nr 
2(T a ) + (V„) = 
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With E b = (T b ) + (V b ) and a = (T a ) + (V a ) one can find 
the kinetic and potential parts explicitly 
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Introducing the pressure p by utilizing the Hugenholtz- 
Van Hove theorem 0] 



X» = E b+P V = E b -v(^) =E b - r f(^) ,7, 
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and eliminating in 1(5} the derivative of E b with the help 
of 0, we obtain (Q]) for the bulk parts. 

With the help of the virial theorem ((31), the electro- 
static potential step at the surface ((3]) can now be ex- 
pressed in terms of the kinetic and potential bulk energy 
and the surface tension as 
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Analogously to [l[ we now employ again the Hugenholtz- 
Van Hove theorem ([7]) and using ([5]) and ([!]) we obtain 



Nn = E b + NA<I>. 



(9) 



This formula is a central result of this manuscript and 
is not found in the literature to the authors' knowledge. 
It describes the chemical potential in terms of the bulk 
energy and the potential step at the surface. The poten- 
tial step at the surface is just the difference between the 
chemical potential and the total energy per particle. 

For infinite (homogeneous) electron systems we can 
shift the surface to infinity. Thus we can now interpret 
the part beyond the bulk energy in the chemical potential 
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Nfi-E b = -(T b ) + -(V b )=pV (10) 



as counting a surface potential step infinitely far away. 
Eq. (fl0|) is exact. The identification of the surface 
potential step with the pressure (|10j) comes from the 
Hugenholtz-Van Hove theorem again. It is physically 
clear that the surface potential step also accounts for the 
pressure. 

Our interpretation also agrees with the expectation 
that the chemical potential should describe the energy 
necessary to add or remove a particle in the system. The 
terms beyond E b expressed in ([I]) can be consequently 
viewed as coming from an imagined surface. Therefore 
adding or removing a particle requires to overcome the 
energy of the bulk as well as the surface potential step. 

This refined picture allows us now to interpret the 
BCS-BEC transition as a competing effect between the 
bulk and surface potential step. For the weak coupling 
(BCS) limit the chemical potential is positive and the 
bulk energy E exceeds the negative surface potential step 
A$. In the strong coupling limit (BEC) the surface part 
can be considered as dominant. The system is unstable 
to the creation of bound states which can Bose condense. 

With this brief summary and reinterpretation of the 
different parts of the chemical potential in the three- 
dimensional electron liquid, we turn now to an integrable 
one-dimensional system. This will allow us to discuss the 
above parts of the chemical potential explicitly. We will 
pay special attention to the different kinetic and poten- 
tial parts as well as the surface potential step. 

In recent work, Wadati and Iida [10( have considered 
the crossover from a BCS-like state to a BEC state in 
a one-dimensional integrable model. The Hamiltonian 
assumed for N particles in a periodic box of length L 
takes the form 
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with m being the mass and c the coupling constant. Fol- 
lowing [Tcj we choose units h = 1 and 2m — 1 for simplic- 
ity. It is useful to define the number 'density' D = N/L 
of particles and the dimensionless coupling parameter 

7=—. (12) 

Wadati and Iida 11011 solve the so-called Gaudin integral 
equation, see also [11| . in the weak coupling BCS regime 
corresponding to I7I <C 1, and hence obtain the ground 
state energy as 
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which according to 
tial parts as 
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FIG. 1: The chemical potential as a function of the total 
energy per particle for the weak (dotted line) and strong cou- 
pling (solid line) limit. 



We now utilize the thermodynamic relation /x = dE/dN 
to get the chemical potential 
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which agrees with the result of [lfj when a misprint is 
corrected in their Eq. (3.25). From (fT5|) with the help of 
(fT3")) we obtain the relation 
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In the limit 7 -> 0, it follows from |Q]| that E B cs/N -> 
tt 2 £> 2 /12, and from HU) that ^ B cs tt 2 D 2 /A. Remov- 
ing 7 in (fl"3"|) and (fl"6|) yields a quadratic equation for the 
chemical potential in terms of the energy, namely 
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where the physical solution is plotted in figure [TJ We see 
that the weak coupling limit 7^2 describes the range 
where fi > and the transition from negative energy to 
positive energy. The latter transition corresponds to the 
transition of bound to free matter. The interpretation 
is that for E < the electrons are bound representing 
Bosons which can form a Bose condensation for fi < 
if the density is high enough in the ground state. In 
order to describe this strong coupling regime, we have to 
consider the other limiting result. 

The scaled ground state energy in the strong coupling 
regime reads 
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which again from® is composed of the kinetic and po- 
tential parts 
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From this we obtain the chemical potential observing 
that D = N/L and 7 = cL/N 
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which corrects some misprint in the last term of Eq. 
(3.8) in [l(|. The curves ^(7) and £(,(7) fit smoothly 
at 7 = —2. We see also that this corresponds nearly 
to the point where the chemical potential changes sign 
which identifies with the BCS-BEC transition point, ex- 
actly being 7^ cs w -1.98 and 7s c trong w -2.07. 

In figure[5]we compare the kinetic and potential contri- 
bution to the chemical potential according to {TJ) calcu- 
lated in ((TH]) and (fT4")l . We see that though the resulting 
chemical potential is smooth at the transition point 7 w 2 
the kinetic and potential energies possess different slopes 
for the BEC and BCS sides. This points towards the 
different correlation mechanisms in both regimes. 

Next we investigate the surface potential step due to 
the new formula (|10[) . We obtain the weak and strong- 
coupling limits 
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which are plotted in figure [3J The surface potential step 
shows a sharp change in the slope at the BCS-BEC tran- 
sition point. While for the strong-coupling limit the sur- 
face potential step is only growing slowly with the energy 
it rises sharply in the weak coupling limit. The diver- 
gence of the strong-coupling limit for 7 > — 2 is present 
in the above formulas for the energy and chemical poten- 
tial as well and is attributed to the approximation (| 18() 
used in (fTU)l where 7 > — 2 exceeds the strong-coupling 
regime. 

The chemical potential /1 in many-body systems con- 
tains fingerprints of phase changes. For the three- 
dimensional correlated homogeneous electron liquid, /1 is 
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FIG. 2: The chemical potential and its kinetic and potential 
parts as a function of the coupling parameter for the weak 
(dotted line) and strong-coupling (solid line) limit. 
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FIG. 3: The surface potential step as a function of the total 
energy per particle for the weak (dotted line) and strong- 
coupling (solid line) limit. The formula (|18p of 10] is not 
applicable for 7 > —2. 



given exactly in terms of the average kinetic and potential 
energies for electrons in eqn. ([1]). The high-density limit 
r s — > is dominated by the kinetic energy. However, as 
the density n is lowered, corresponding to increasing r s , fi 
passes through zero, which is related to surface phenom- 
ena. Utilizing the Budd-Vannimenus theorem we have 
been able to identify the difference between the chemical 
potential and the total energy per particle as the surface 



potential step assuming that we have an infinite flat sur- 
face of half space. This allows us to view the sign change 
of the chemical potential as a competing effect of total 
energy and surface potential step. 

Though we have restricted our discussion to a uniform 
electron liquid, for /i sufficiently negative the liquid is 
unstable against Wigner electron crystallization [12j | . 

Turning finally to the one-dimensional model of Wa- 
dati and Iida [lfj possessing a BCS-BEC crossover, it 
has been demonstrated that the chemical potential fj, can 
be expressed solely in terms of the ground-state energy 
per particle. The coupling constant of the electrons in 
the many body Hamiltonian no longer appears explicitly. 
We calculated the surface potential step according to our 
finding (flQ|) and observed that the BCS-BEC transition 
is characterized by a change in the slope of the surface 
potential step with the ground state energy per particle. 
This could be viewed as an alternative measure for the 
BCS-BEC transition. 

It would, of course, be interesting if an analytically 
tractable model in higher dimensions could be devised in 
which fi = n(E/N) could be obtained, showing a BCS- 
BEC crossover. 
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